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We show numerically that the ‘deconfined’ quantum critical point between the Neel antiferro- 
magnet and the columnar valence-bond-solid, for a square lattice of spin-l/2s, has an emer¬ 
gent SO(5) symmetry. This symmetry allows the Neel vector and the valence-bond-solid order 
parameter to be rotated into each other. It is a remarkable 2-|-1-dimensional analogue of the 
SO(4) = [SU(2) X SU(2)]/Z2 symmetry that appears in the scaling limit for the spin-1/2 Heisen¬ 
berg chain. The emergent SO(5) is strong evidence that the phase transition in the 2-l-lD system 
is truly continuous, despite the violations of finite-size scaling observed previously in this problem. 

It also implies surprising relations between correlation functions at the transition. The symmetry 
enhancement is expected to apply generally to the critical two-component Abelian Higgs model (non¬ 
compact CP'^ model). The result indicates that in three dimensions there is an SO(5)-symmetric 
conformal field theory which has no relevant singlet operators, so is radically different to conventional 
Wilson-Fisher-type conformal field theories. 


Many condensed matter systems show higher symmetry 
in the infrared than they do in the ultraviolet. The liquid- 
gas critical point is a classical example: although there is 
no microscopic Z2 symmetry exchanging liquid-like and 
gas-like configurations, the fixed point has an emergent 
Z2 symmetry and is simply the Ising fixed point. Mi¬ 
croscopically this fixed point is perturbed by operators 
which break the Z2 symmetry, but it nevertheless gov¬ 
erns the critical behaviour because these perturbations 
are irrelevant under the renormalisation group. 

To reach this critical point two variables, say tem¬ 
perature and the pressure, must be tuned. The spin- 
1/2 Heisenberg chain provides an example of emergent 
symmetry without such fine-tuning in a quantum set¬ 
ting. The ground state of this model is well known to be 
critical. Its microscopic symmetries are SU(2) spin ro¬ 
tations, together with spatial symmetries. However the 
scaling limit of the spin-1/2 chain is the SU(2)i Wess- 
Zumino-Witten conformal field theory [T], and this has 
an SO(4) = [SU(2) x SU(2)]/Z2 symmetry which is much 
larger than the global symmetry present microscopically. 

Physically, this arises as follows [5]. The Neel vector N 
has three components. There is also a spin-Peierls order 
parameter ip which distinguishes between the two differ¬ 
ent patterns of dimer (singlet) order, and which changes 
sign under appropriate reflections or translations. We 
may form the 4-component superspin $ = (N,ip), and 
the emergent SO (4) corresponds to rotations of this vec¬ 
tor. Although the dimer and Neel order parameters are 
utterly inequivalent microscopically, a symmetry between 
them arises in the infra-red. Technically, this again relies 
on the SO(4)-breaking perturbations of the conformal 
field theory being irrelevant or marginally irrelevant. 

Naively one might expect this phenomenon to be spe¬ 
cial to one spatial dimension, where the enlarged sym¬ 
metry is related to special properties of 2D conformal 
invariance (the doubling of conserved currents [1]). We 
show here however that an analogous symmetry enhance¬ 


ment occurs for the spin-1/2 magnet on the square lat¬ 
tice, at the celebrated ‘deconfined’ quantum critical point 
[3H5]. This is a transition between the antiferromagnet- 
ically ordered Neel state and a columnar valence-bond- 
solid (VBS), and is reached by tuning a single parameter. 
For example, the nearest-neighbour Heisenberg model 
can be driven into the VBS using either a four-spin in¬ 
teraction [7] or a next-nearest-neighbour exchange mm- 
The emergent symmetry we put forward is an SO(5) 
which mixes the components of the Neel vector N, which 
has three components, and the VBS order parameter (p, 
which has two. We test it by examining the joint proba¬ 
bility distribution for these quantities. 

Numerically, the critical behaviour can be studied effi¬ 
ciently with a 3D classical loop model [9] , and we use this 
approach here. The order of the transition has been con¬ 
troversial as a result of violations of conventional finite- 
size scaling [lOHIK] which we discussed in detail previ¬ 
ously [H]. We will return to this below, arguing that the 
present results support the continuity of the transition. 

The Neel-VBS transition is usually described with the 
non-compact CP^ (NCCP^) Lagrangian [4l[T6]. 

£ = |(V - iA)z\^ A- k(V X Af + pl\z\^ + (1) 

The two-component bosonic spinon field z = (21,02) is 
related to the Neel vector iV by iV = z"^az. The U(l) 
gauge field A^ is related by duality to the VBS order 
parameter 0 = {ipx,ipy) which distinguishes the different 
columnar singlet patterns mm Although we will use 
the language of the Neel-VBS transition, our conclusions 
apply more generally to the above field theory, and in¬ 
dicate that it flows to an SO(5)-symmetric fixed point 
at the critical value of p.. In the language of this 3D 
gauge theory, the VBS order parameter is the operator 
M. = ipx + iipy which inserts a Dirac monopole in 
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so(5) symmetry cannot be made explicit in the for¬ 
mulation of Eq. Fortunately, Senthil and Fisher |20) . 
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FIG. 1. The joint probability distribution P{Nx,^x), after 
rescaling N^; and ifix, to have unit variance, in a critical system 
of size L = 100. Upper plane shows contour plot. 

building on work of Tanaka and Hu [21) . have argued 
that an alternative field theory describes the Neel-VBS 
transition and is equivalent to Eq. This is a nonlinear 
cr-model (NLcrM) for the five-component superspin 

^ = {N^,Ny,N^,^^,(Py), ( 2 ) 

augmented with: (i) anisotropies that break the global 
symmetry from SO(5) down to the spin rotation and spa¬ 
tial symmetries present microscopically; and (ii) a topo¬ 
logical Wess-Zumino-Witten term at level one [5^ |23] , 
which is analogous to that in the GET for the spin chain. 
The leading anisotropy plays the role of the mass term 
in Eq. it drives the transition between the Neel and 
VBS ordered phases. 

The NLcrM formulation makes the emergent symme¬ 
try a more natural possibility, since it could arise at 
the critical point if all the higher anisotropies happen to 
be renormalisation-group irrelevant. We will discuss be¬ 
low the phase diagram for the NLcrM (with WZW term) 
which is implied by this conjecture. 

In previous work we characterised various observables 
at the deconfined transition in detail, using a three- 
dimensional loop representation to reach system sizes up 
to L = 512. See Ref. [3] for details of the model, which 
is in the Neel phase for coupling J < Jc and the VBS 
phase for J > Jc, with Jc = 0.088501(3). We found a 
remarkable similarity between the critical Neel and VBS 
correlation functions. The anomalous dimensions deter¬ 
mined from the correlators at separations r L are 
VNeei = 0.259(6) and ryvBS = 0.25(3) [23]; the two corre¬ 
lators also behave similarly in the regime r ^ L (despite 
the scaling violations discussed in Ref. 0 ). This suggests 
searching for an emergent SO(5) symmetry that would 
explain these apparent coincidences. 

Probability distribution. Consider the joint distribu¬ 
tion for the Neel and VBS order parameters in a system 
of linear size L. If SO(5) emerges, then this will be a 


function only of -|- at the critical point (af¬ 

ter a trivial rescaling of ip). Spin rotation symmetry of 
course already guarantees that the distribution depends 
on N only via V^. Also, while microscopic spatial sym¬ 
metry only allows (p to rotated by multiples of 7r/2, it is 
well established numerically that symmetry under con¬ 
tinuous U(l) rotations of tp emerges close to the transi¬ 
tion |25l [26]. This was checked for the present model in 
Ref. |9] (see also App. and see Ref. m for related 
phenomena). The crucial point is therefore whether the 
distribution is invariant under U(l) rotations that mix a 
component of p with a component of N. 

Let the standard deviations of and px be denoted 
CTjv and Oy, respectively, and use a tilde to denote quan¬ 
tities rescaled to have unit variance: Nx = Nx/ctn and 
Px = Pxla^- Fig. 0 shows the joint probability distribu¬ 
tion for these quantities at the critical point J = J^ va. a 
system of size L = 100. The visual evidence for emergent 
symmetry between Neel and VBS is striking. 

Before turning to a quantitative analysis of the distri¬ 
bution, a basic test is that the variances ajq and of 
the order parameters depend on system size in the same 
way at criticality |2H]: be. that is L-independent 

at Jc- In Fig. [^this is confirmed to high precision over 
a wide range of lengthscales. Plots of a^pjat^ versus J 
for different L-values cross at Jc [23] (the value of 
at Jc depends on the nonuniversal normalisation of the 
lattice operators). 

For a quantitative analysis of the probability distribu¬ 
tion we examine the moments 

F/= (r“cos(£6l)), (3) 

where {Nx,Px) = r (cos 0, sin 0). Emergent symmetry 
requires these to vanish for ^ > 0. We have computed F 2 



FIG. 2. Main panel: variance ratio (Jp/aN plotted against 
J for various L. Gurves cross at Jc as expected from SO(5) 
symmetry. Inset: same quantity as a function of L for several 
J around Jc ^ 0.0885 (key in Fig. j^. 
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and for large sizes: 

F^ = (N^-ipt), F^ = (^N^-6Niipl + ^t). (4) 

Fig.0 shows these as a function of L at and close to the 
critical point. Both are consistent with zero for J = 
over the whole range of L. The expected values in the 
adjacent phases (including the regime of weak VBS order, 
where there is an effective U(l) symmetry for rotations 
of If) are also indicated in the figure (details in App. H- 

In Fig. 1^ we show F 2 as a function of J: note the 
very clearly defined crossing at J = Jc, F 2 = 0- Further 
moments are shown for L = 100 in App. It should 
be noted that the critical distribution is markedly non- 
Gaussian, with nonvanishing higher cumulants such as 
[{N^} - 3(iV2)2]/(iVj)2 = -0.7549(13) (for L = 100). 

Equalities between scaling dimensions. In addition to 
the equivalence between Neel and VBS vectors (mani¬ 
fested in the joint distribution and the anomalous dimen¬ 
sions), SO (5) has consequences for operators transform¬ 
ing in higher representations. Take the leading operators 
in the symmetric two- and four-index representations: 

OilL = " Catcd- 

The subtractions |30| ensure the operators transform irre- 
ducibly. is relevant, with scaling dimension X 2 < 3. 
In fact a component of is the operator Oj which 
drives us through the Neel-VBS transition as we vary J, 
by favouring one or the other order {Oj therefore plays 
the role of the mass term in Eq. : 

= ( 5 ) 

Remarkably, various a priori unrelated operators share 
the same scaling dimension X 2 since they are also com- 
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FIG. 3. Moments F 2 (main panel) and F* (inset) defined in 
Eq.H shown as a function of L for a few values of J at and 
close to Jc ~ 0.0885. Dashed lines indicate expected values 
in the Neel phase, at the SO(5)-symmetric critical point, and 
in the U(l) and Z 4 -symmetric regimes of the VBS phase. 



FIG. 4. Fi (Eq. § as a function of J for various sizes, showing 
a crossing at F 2 = 0, J = Jc as expected for SO(5) symmetry. 

ponents of These include the spin-quadrupole mo¬ 
ments, NaNb — dat,iV^/3, and the relevant [26] opera¬ 
tors fa^b — which in the NCCP^ language in¬ 

sert ‘strength-two’ monopoles m- Even more oddly, the 
same scaling dimension controls ipaNb, despite the fact 
that microscopically this operator is as dissimilar from 
Oj as possible — fa^b transforms under both spin and 
spatial symmetries, while Oj is invariant under them. 

To test these predictions. Fig. shows data for the 
two-point functions of Oj, ipxNz, and ipxfy, or rather 
lattice versions of these operators. (See App. |^for defi¬ 
nitions and Ref. jS] for a general discussion of correlation 
functions at the critical point.) Note the striking simi¬ 
larity of the three curves, as expected from SO(5). The 
slopes at r ^ 10 are around xf^ ~ 1.5, but this effective 
exponent may be strongly affected by finite size effects. 
(Various scaling dimensions, including the two-monopole 
dimension |31j . are known at large n in the SU(n) gen¬ 
eralisation of Eq. , and show that a symmetry 

between Neel and VBS cannot persist in this limit.) 

allows us to write both a subleading opera¬ 
tor which breaks the symmetry between Neel and VBS 
(X]a=i X]fc =4 which breaks the remaining 

symmetry for ip down to fourfold rotations (X]a =4 ^aaaa)- 
Therefore it is possible that the same irrelevant exponent 
controls finite-size corrections to both types of symmetry 
enhancement (see also App. H- 

Nonlinear a-model. The NLctM description of the de- 
confined critical point proposed in Ref. [20] is 

5,, = / Q (V$)' + ^ + 5wzw (6) 

where 5wzw is a topological Wess-Zumino-Witten term 
at level one (associated with the homotopy group 
7r4(S^) = Z of the target space). Physically, this term 
ensures that a vortex in the VBS order has an unpaired 
spin-1/2 at its core [TT] HHI- The Oi are the various 
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anisotropies, some discussed above, that break SO(5) 
symmetry down to the microscopic physical symmetry. 

Suppose that the Neel-VBS transition is continuous, 
with emergent SO(5) symmetry, and that the NLctM is 
a viable description of this transition. Since the criti¬ 
cal point is reached by tuning a single parameter, there 
is only a single RG-relevant coupling in Eq. namely 
the anistotropy Oj of Eq. Therefore the SO(5)— 
symmetric NLcrM, with no anisotropies, has a nontrivial 
infra-red stable fixed point controlling a power-law corre¬ 
lated phase: see Fig. This fixed point then also governs 
the deconfined transition. (The NLctM also has a stable 
ordered phase even in the presence of the WZW term.) 

The phase diagram of Fig. is counterintuitive, as we 
are used to 3D cr-models with an ordered phase, a triv¬ 
ial disordered phase, and an unstable critical point in 
between. But that conventional picture applies to the 
NLctM without a WZW term. The present results in¬ 
dicate that the WZW term causes the trivial disordered 
phase to be replaced with a stable critical phase. While 
the absence of a trivial disordered phase may seem sur¬ 
prising, it is in fact a necessary feature of any field theory 
which faithfully describes the spin-1/2 magnet, as argued 
in Ref. [20] . By the 2D generalisation of the Lieb-Schultz- 
Mattis theorem [Ml [37], the magnet cannot be in a triv¬ 
ially disordered phase with a unique, symmetric, gapped 
ground state (it can be critical, or topologically ordered, 
or spontaneously break a symmetry). Thus the field the¬ 
ory should not have such a trivial phase either |41j . This 
argument does not by itself imply the existence of a sta¬ 
ble critical phase, only the absence of a trivial fully dis¬ 
ordered one. (See Refs. |5DJ [35] for related discussions of 
the SO(4) case.) 

The emergent SO(5) therefore indicates the existence 
of a 3D SO(5)-symmetric GET which is radically unlike 
standard Wilson-Fisher CFTs, in that there are no rel¬ 
evant singlet operators. It would be very interesting to 
investigate this using the conformal bootstrap [42] , mak¬ 
ing use of numerical estimates for operator dimensions 



FIG. 5. Correlators Go{r) for operators ifxNz, ipxg>y and Oj 
in a system of size L = 100. (The Go have been normalised 
to agree at r = 10.) Inset: Goj for various L. 
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FIG. 6. Conjectured phase diagram for NLctM with WZW 
term in presence of full SO (5) symmetry. The fixed point on 
the left also governs the deconfined critical point, where SO(5) 
is emergent. Moving away from Jc introduces the relevant 
symmetry-breaking perturbation Oj, leading to the Neel or 
VBS phase. 


[9]. This should be simpler [43] than studying the criti¬ 
cal NCCP^ model without assuming SO(5). 

Scaling violations. The deconfined critical point shows 
strong violations of finite-size scaling. We argued in 
Ref. [9] that these are not simply large corrections to 
scaling of the conventional type (i.e. from irrelevant 
or marginally irrelevant operators), but should instead 
be attributed either to a first order transition with an 
anomalously large correlation length or to a genuine crit¬ 
ical point with unconventional finite-size scaling due to 
a dangerously irrelevant variable. The present results 
strongly suggest the second scenario — a genuine con¬ 
tinuous transition. This is because a critical point is the 
only natural explanation for the emergent SO(5) symme¬ 
try, which we have tested here to high precision. This 
symmetry therefore provides long-sought direct evidence 
that the deconfined transition is continuous. 

In more detail: if one attempts to account for the data 
in terms of an anomalously weak first order transition 
(i.e. without postulating a genuine 3D critical point), 
one is led to a scenario where the apparent critical be¬ 
haviour is due to a ‘nearby’ fixed point at spacetime di¬ 
mension slightly below three |9]. While this scenario can 
potentially explain pseudo-critical behaviour up to an ex¬ 
tremely large lengthscale, it cannot account for the emer¬ 
gent SO(5) symmetry — this makes sense only for a 3D 
fixed point. While we can consider the NCCP^ model in 
arbitrary dimension, the operator (p (interpreted for ex¬ 
ample as a monopole insertion operator) is special to 3D, 
and is required for construction of the SO(5) superspin. 

Assuming therefore that the transition is continuous, 
a possible explanation for the scaling violations is that a 
dangerously-irrelevant variable is required to cut off the 
fluctuations of a zero mode of the field [9]. (This would be 
analogous to theory above four dimensions, where the 
quartic term is irrelevant, but cannot be set to zero since 
it is needed to prevent divergent fluctuations of i/’s zero 
mode [44].) This may suggest that an alternative field 
theory description exists which is more natural than the 

NLctM [33] ■ 

Future directions. It would be interesting to investi¬ 
gate consequences of the emergent symmetry for finite 
temperature behaviour, as well as to look for signs of it 
using methods complimentary to Monte Carlo such as 
exact diagonalisation or DMRG. The present results also 
motivate further analysis of SO(5)-symmetric 3D CFTs 
— for example with the conformal bootstrap — and a 
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more general investigation of the role played by WZW 
terms in field theories above two dimensions. For exam¬ 
ple, is there an analogous SO(6)-symmetric CFT in 4D? 
Finally, we note that the critical behaviour at the decon- 
fined transition remains perplexing, and deserves further 
investigation. 
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Appendix A: Lattice definitions of operators 

We refer to the lattice field theory of Ref. [5], de¬ 
scribed there for J = 0 but modifiable for J ^ 0. The 
Neel vector N is defined on each link hy N = z^az 
with |zp = 2 (these normalisations are purely a mat¬ 
ter of convenience). A graphical expansion maps the 
lattice field theory to a partition function for loops tak¬ 
ing two colours, red and blue [33]. Inserting operators 
on the links modifies the graphical expansion. Using 
= (2/3)5cd(iJci — ^c 2 ), where Tr is the inte¬ 
gral over zs and Tt ZcZ’^ = Scd, one may check that is 
simply the operator which measures the colour of a link 
in the loop gas: = (2/3)x, where x is ±1 depending 

on whether the link is red or blue. Therefore we obtain 
the probability distribution for the z-component of the 
uniform magnetisation, Ad°* by measuring the number of 
red links in the Monte Garlo simulation. We may mea¬ 
sure ip simultaneously, allowing construction of the joint 
probability distribution. This is what we have done for 
system sizes L < 100. (In the main text the component 
we singled out was labelled rather than Nz, but of 
course by symmetry there is no difference.) 

For larger system sizes we do not have data for the full 
probability distribution of link colours. However we have 
data for low moments of the loop length distribution, 
and this is sufficient to construct the moments and 
shown above. We define, in a given configuration, 

Sk = ^ (length of loop)*. (Al) 

loops 

The necessary relations follow straightforwardly from 
= (2/3)EunksX- We express ((fV‘°‘)^) as a 
correlation function of x operators in the loop gas. Sum¬ 
ming over the possible colourings of a given loop config¬ 
uration gives a nonzero result only if there are an even 
number of x insertions on each loop. Taking into account 
the possible ways of assigning ys to loops, we have 

((fVfT) = (2/3)^(52), (A2) 

((Af‘)4) = (2/3)^ {3Si - 2 S 4 ). (A3) 


pa 

a = 0 

a — 2 

a = 4 

a — & 

1 = 2 

0.00035(17) 

0.000(4) 

0.00(2) 

0.01(7) 

1 = 4 

0.0002(7) 

-0.001(2) 

-0.004(8) 

-0.02(3) 

1 = 6 

0.0002(5) 

0.000(1) 

-0.003(5) 

-0.02(2) 


TABLE I. Moments in the {Nx,ipx) plane for L = 100 
(taken at J = 0.0885). 



3 -3 


FIG. 7. The joint probability distribution P(ipx,Px), after 
rescaling and Py to have unit variance, in a critical system 
of size L = 100. Upper plane shows contour plot. 

Importantly, these relations remain true if powers of ip 
are inserted on the left and right hand sides. This allows 
us to construct the abovementioned moments. 

Note that in order to give optimal statistics, the uni¬ 
form magnetisation/VBS order is always defined by an 
integral over the three-dimensional sample — i.e. space- 
time, in the quantum language — rather than a two- 
dimensional ‘spatial’ slice. 

The 3D ‘L’ lattice [46] on which the model is defined 
is bipartite, with defined at nodes of the A sublat¬ 
tice and py at nodes of the B sublattice, and the Boltz¬ 
mann weight couples nearest neighbour p components on 
the same sublattice |9] . The operators whose correlation 
functions are plotted in Fig. are all defined at a node, 
and the correlation functions are evaluated for pairs of 
nodes on the same sublattice, separated parallel to a co¬ 
ordinate direction. We define Oj{r) as the sum of the 
eight terms in the lattice energy which involve the node 
r, and subtract the average so that {Oj) = 0. The op¬ 
erator [pxPy]{r) is defined as the sum of pa;(r)p{r') over 
the four nearest neighbours r' of r (see Fig. 1 in Ref. |9|) 
and [pz:Nz\{r) is defined as the sum of Px{r)Nz{l) over 
the four links I connected to r, with Nz defined as above 
in terms of link colours. Each of these operators is there¬ 
fore supported on star: we take the separation of the two 
points in the correlator to be perpendicular to the plane 
of the stars. 
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Appendix B: Additional data for probability 
distribution 

Table |T] shows data for a variety of the moments 
defined in Eq. [^in a system of size L — 100 at J = 0.0885. 
These quantify the U(l) symmetry of the distribution in 
the {Nx,ipx) plane. All are consistent with zero. 


L 

a = 0 

a = 2 

a = 4 

a = 6 

32 

-0.0086(9) 

-0.035(3) 

-0.166(10) 

-0.88(5) 

64 

-0.005(1) 

-0.021(3) 

-0.095(11) 

-0.48(5) 

100 

-0.0028(14) 

-0.014(4) 

-0.069(15) 

-0.35(7) 


TABLE II. Moments in the {ipx,^y) plane, quantifying 
fourfold anisotropy for the VBS, for L = 32,64,100. (Note 
that moments grow with a for trivial reasons.) 

We also examine the symmetry properties of the prob¬ 


ability distribution in the (ipxT^y) plane. This is shown 
in Fig. for a critical system of size 100. Again it is 
U(l) symmetric to an excellent approximation. Interest¬ 
ingly though, the finite-size corrections to this U(l) are 
larger than those for the U(l) in the (ATj,, (fx) plane (but 
still small). Table |n] quantifies fourfold anisotropy for 
the VBS at the critical point via the moments (de¬ 
fined analogously to but in the {Ipx, ^y) plane). They 
decrease with system size roughly as L~‘^ with c ~ 0 . 8 . 

In Fig. of the main text the straight lines indicate 
the expected values of the moments ^ad within 
the phases. These values follow straightforwardly. The 
vector which is in disordered phase has a Gaussian distri¬ 
bution, so its moments follow from Wick’s theorem. The 
other order parameter can be treated as fixed in magni¬ 
tude and averaged over symmetry-equivalent directions. 
In the Neel phase this gives F^ = F^ = —1.2. Our 
choice of coordinates is such that deep in the VBS phase 
If oc (±1,±1), which gives F 2 = 2, F 4 = —2. For weak 
VBS order there is a large range of L where U(l) sym¬ 
metry for If survives, and F 2 = 1.5, F^ = —1.5. 
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